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Abstract 

We present an ah initio study of the recently discovered superconductivity 
of boron doped diamond within the framework of a phonon-mediated pairing 
mechanism. The role of the dopant, in substitutional position, is unconven- 
tional in that half of the coupling parameter A originates in strongly localized 
defect-related vibrational modes, yielding a very peaked Eliashberg a^F{uj) 
function. The electron-phonon coupling potential is found to be extremely 
large and Tc is limited by the low value of the density of states at the Fermi 
level. 
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EFFET ISOTOPIC DU BORE: SUGGERER EXPERIENCE !! 

In a recent paper^, the superconductivity of boron doped (B-doped) diamond was ev- 
idenced experimentally with a transition temperature Tc of 4 K. This work follows the 
discovery of superconductivity in doped silicon clathrates ^'^ (Tc=8 K), a cage-like sihcon 
material which crystallizes in the same sp^ environment as the diamond phase. Further, 
the superconductivity of carbon-based clathrates was predicted on the basis of ab initio 
simulations^'^. 

Even though the reported temperatures are rather low, the superconducting transition 
of column IV semiconductors is of much interest as it concerns very common materials. 
Further, in the case of the carbon clathrates, it has been predicted that the electron-phonon 
potential Vep=A/N(Ei;') (where N(Ei^) is the density of states at the Fermi level), may be 
extremely large'^''^. Finally, the superconductivity of hole-doped strongly bonded covalent 
systems is reminiscent of that of the MgB2 compound ^. 

As compared to the cage-hke clathrate phase, diamond is difficult to dope as the network 
is much denser. This explains why the superconductivity of doped diamond was never 
observed, as only quite recently a doping in the limit of a few percent could be achieved. 
As described in Ref. 1, it is only for a boron concentration c > 10^^ cm~^ that a degenerate 
semiconductor or metallic character can be observed. 

In two recent work^'^, a first theoretical study of the superconductivity of B-doped dia- 
mond within an elegant virtual crystal approximation was proposed. In this approach, the 
crystal is made of "averaged carbon-bore atoms" (nuclear charge Z=(l-x)Zc + xZ^) within 
the standard 2 atoms/cell structure of diamond, and B atoms are not considered explicitly. 
As such, the superconductivity of B-doped diamond was analyzed in terms of Bloch-like 
delocalized electron and phonon "average" states. 

In this work, we study by means of ab initio calculations the electron-phonon coupling in 
B-doped diamond. We consider explicitly a large supercell with one B atom in substitution. 
The metallic character of the doped system is evidenced and the deviations from a simple 
rigid band model are discussed. We show further that the boron atoms yield very localized 
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vibrational modes which contribute to half of the electron-phonon coupling parameter A, in 
large contrast with the picture conveyed by the virtual crystal approximation. The value of 
A is calculated ab initio, yielding an unusually large electron-phonon coupling potential Vgp 
or deformation potential D. 

Our calculations are performed within the local density approximation to density func- 
tional theory® (DFT) and a pseudo-potential^ plane-wave approach. A 55 Ryd energy cut-off 
and a (4x4x4) Monkhorst-Pack^° sampling of the Brillouin zone (BZ) showed good conver- 
gence for structural relaxations. The phonon modes and the electron-phonon (e-ph) coupling 
matrix elements were obtained within the framework of perturbative DFT^^. Due to the 
computational cost, phonons were calculated at the q = F point only. This is equivalent to 
a (3x3x3) sampling of the undoped diamond BZ. 

For the calculation of the average e-ph matrix elements gf^^ over the Fermi surface, 
required to evaluate A: 

A = N{Ef)Vep = 2N{Ef) J2 < > /hcug, (1) 

< >=fi^T. (2) 

a much larger (8x8x8) k-point sampling was used. In Eqn. 2, the origin of the electronic 
energies (efe„) is taken to be the Fermi energy Ep. 

Our system consists of a diamond (3x3x3) supercell (54 atoms) with one B atom in 
substitution (noted BC53 in what follows). We are thus in a 1.85% doping limit, close to 
the 2.5 ±0.5% studied experimentally. Upon structural relaxation, the B-C bond length is 
1.57 A as compared to the 1.54 A for the C-C bonds in undoped diamond (LDA vahics). 
The distance between second and third carbon neighbors is reduced to 1.52 A in order to 
accommodate the larger B-C bonds. 

We first explore the electronic properties and compare the band structure of BC53 with 
that of diamond (see Fig. 1). As expected, we are clearly in a metallic regime with the Fermi 



level located at ~ 0.54 eV below the top of the valence bands. Even though Fig. 1 seems 
to confirm the model of rigid band doping around E^, clear differences exist in particular 
at the zone boundaries where degeneracies are left due to the reduction of the periodicity 
under doping. We note that in most directions, the band structure of pure and B-doped 
diamond start to significantly depart from each other precisely below E^^^^. 

We also calculate the total density of electronic states (eDOS) projected onto the atomic 
orbitals. Even though the eDOS projected onto the B orbitals (cDOSb) is small, it is 
instructive to plot (Fig. lb) the ratio NaXeDOSs/eDOS, where Na=54 is the number 
of atoms per unit cell. Clearly, close to Ep, the relative weight of B orbitals into the 
wavef unctions is much larger than the weight on any given carbon atom. As expected, this 
ratio increases where the bands of the ideal and doped diamond differ most significantly. 
An analysis of the projection of the eigenstates onto the atomic orbital basis shows that the 
leading coefficients are the {px,Py,Pz) B orbitals^^. Such results indicate that even though 
degenerate with the diamond Bloch states, the wavefunctions of interest for e-ph coupling 
exhibit some degree of localization around the B atom. 

We now study the phonon states at F. The related DOS (pDOS) is shown in Fig. 2b and 
compared in Fig. 2a to the one of undoped diamond As compared to the virtual-crystal 
approximation results, we do not observe the very large softening of the zone-center phonon 
modes predicted in Ref. 6. However, the splitting of the highest diamond mode (Fig. 2a) 
clearly results in transferring spectral weight to lower energy. This can be interpreted on the 
average as a softening of the optical frequencies. In particular, zone-center modes emerge 
around 1200 cm~l (Fig. 2b) which are consistent with the appearance of new broad peaks 
in Raman experiments around 1220 cm~^^'^^. 

More insight can be gained by calculating the pDOS projected onto the B atom and its 
4 C neighbors (thick hne). The most salient feature is that the new modes (3-fold) at 1024 
cm~^ (vertical dotted hues in Fig. 2) project to more than 50% on these 5 atoms. The 
analysis of the associated eigenvectors clearly indicates that these three modes correspond 
to the stretching of the B to neighboring C atom bonds. The displacement of the B atom is 
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~ 22% larger than that of C neighbors as expected from the difference in masses. Therefore, 
the dopant atoms induce vibrational modes which are very localized in space. This is a 
significant deviation from a virtual crystal approach where only extended phonon states can 
be evidenced. The influence of such modes on the e-ph coupling is now shown to be crucial. 

We follow equations 1-3 to calculate A. With our hmited q-point samphng^^, we find 
that X{q = r) is rather sensitive to the gaussian broadening used to mimic the 5-functions 
in Eq. 2. However, such a variation is primarily related to the variation of the nesting factor 
n(q) (that can be obtained by setting gqi,{knn') to unity in Eq. 2). This factor n(q) is a 
measure of phase space availability for electron scattering by phonons. Indeed, one verifies 
that the quantity \{q)/n{q) (q=r) is much more stable and displays a plateau over a large 
energy broadening range. We thus define an averaged A by rewriting: 

\ = j d'qX(q)^ I d\n{q)\{q)/n{q) 
- A(r)/n(r) j d\niq) = iV(E^)2A(r)/n(r) 

Such an approximation relies on the usual assumption that the e-ph matrix elements are 
varying much more smoothly as a function of q than the nesting factor. 

We find A=0.43 ±0.01. Such a value of A is much smaller than that of MgB2 (in the ab- 
sence of anharmonic effects^^, A(MgB2) ~ 1). Assuming now that the McMillan expansion-*^^ 
is valid in the limit of degenerate semiconductors: 

-1.04(1 + A) 
A-/x*(l + 0.62A)_ 

and with a logarithmic averaged phonon frequency cuiog of 1023 cm~^, the screening pa- 
rameter //* needed to reproduce the 4 K experimental value of T^ is found to be of the 
order of 0.13-0.14^^. This is a very standard value, indicating that phonon-mediated pairing 
mechanisms can adequately account for the superconducting transition in impurity-doped 
diamond. 

With N(Ef)=20.8 states/spin/Ry/Cell, the e-ph potential Vep is 280 meV ± 10 meV. 
Such a value is comparable to what was found in the Li-doped carbon clathrates (V ep ~ 250 
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meV^) and much larger than the 60-70 meV found for the fuUerides This was analyzed 
in terms of the amount of sp^ character of the C-C bonds^. For sake of comparison with 
Refs. 6, 7, we calculate the deformation potential D such that A=N2D^/Ma;^, where N2 is the 
eDOS in eV/state/spin/"2 atoms cell" . With N2=0.057 and a;=1023 cm-\ one finds D = 27 
eV/A~^, larger than the D ~ 21 eV/A~^ value found in the virtual crystal approximation^'^. 
This value is more than twice as large as the one calculated for MgB2, indicating that it is 
primarily the low density of states at the Fermi level, related to the 3D nature of diamond, 
that is responsible for the low value of A and Tc. 

Better q-point sampling and the influence of anharmonicity^'^^ will certainly change 
the value of A. Beyond the numerical estimates, it is important to look at the spectral 
decomposition of A, that is the Ehashberg function a'^F{uj) (Fig. 2c): 

Q;^F(a;) = ^"^UJqu^uSiLO - LUqu) (4) 

By comparison with the pDOS, it appears that a significant contribution to A comes from 
the B-related phonon modes. As a matter of fact, these 3 modes contribute to about 50% of 
the coupling^^. In particular, the resulting Eliashberg function is very peaked at the impurity 
vibrational modes energy, displaying a S-\ike shape. Therefore, not only the density of states 
N(Ei7'), but also the strength of the e-ph coupling, can be taylored by changing the dopant 
concentration and/or chemical type. This mechanism offers thus an additional flexibility in 
trying to increase T^;. The situation is very different from the one of fuUerides or clathrates 
where the low frequency phonon modes associated with the dopant hardly contribute to A. 

To strengthen that point, we look at the superconductivity of diamond within a rigid 
band modeP^. We plot in Fig. 3 the evolution of A as a function of the position of the 
Fermi level. For Ej located ~ 0.5-0.6 eV below the top of the valence bands, we find a 
A of ~ 0.24, that is roughly 55% of what we find in the BC53 case. Even though very 
qualitative, this computer experiment substantiates the conclusion that B-related modes 
contribute significantly to the e-ph coupling. Further, Fig. 3 clearly suggests that larger A 
could be achieved under larger doping conditions due to the quadratic rise of N(Ei7) as a 
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function oiEp- 

It is certainly accidental that just half of the coupling originates in B-modes as the weight 
of the eigenstates on the dopant should scale as \/Na under varying concentration (neglecting 
the electronic localization effect described above). In addition, the simple structural model 
studied here does not account for the effect of randomness, nor the possibility for B atoms to 
cluster or adopt other configurations. It is clear however that the role of the dopant cannot 
be interpreted solely as providing charges to the network, nor as softening the phonon modes 
by depleting the bonding states at the top of the valence bands. In particular, such results 
reveal that any attempt to explain the superconductivity of doped diamond within a rigid 
band model or the virtual-crystal approximation should be taken with care. The importance 
of the impurity related vibrational modes is expected to be very general and should apply 
to any semiconductor/insulator doped by strongly bound atoms. 

In conclusion, we have shown that the role of boron in doped diamond is unconventional 
in that a significant fraction of the coupling coefficient A originates in vibrational modes 
very localized on the defect centers. With a A of 0.43, much smaller than that of MgB2 ~ 1, 
Tc remains small. Even though the electron-phonon potential is extremely large, the 3D 
nature of the network reduces the density of states at the Fermi level. This invites to study 
the case of doped diamond surfaces where both the contraction of the reconstructed bonds 
and the 2D nature of the surface states may lead to much larger Tc^^. 
Acknowledgements: Calculations have been performed at the French CNRS national 
computer center at IDRIS (Orsay). The authors are indebted to region Rhone-Alpes for 
financial support through the "Thematique prioritaire programme 2003-2005" in Material 
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FIGURES 

FIG. 1. (a) Details of the top of the valence bands along high-symmetry directions of the FCC 
BZ for undoped diamond in the (3x3x3) cell (dotted lines) and the doped BC53 system (full lines). 
The horizontal line is the position of the Fermi level in the doped case. The two band structures 
have been aligned at the top of the valence bands, (b) Ratio of the density of states projected onto 
the B-orbitals to the total density of states per atom in the unit cell (this ratio has been set to 
zero in the band gap). 

FIG. 2. Phonon density of states pDOS for (a) bare diamond with modes on an unshifted 
(3x3x3) grid (b) the boron doped BC53 cell with modes at F. In (c) the spectral function a^F{u) 
is represented for modes at F. Energies arc in cm~"^ and the y-axis magnitude is arbitrary. A 5 
cm~^ broadening has been used. The vertical dotted lines indicate the energy position of the 3-fold 
B-related vibrational modes. 

FIG. 3. Evolution of A for diamond, doped within a rigid band model (see text), as a function 
of the Fermi energy position E^r. The energy reference has been set to the top of the valence bands. 
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We present an ab initio study of the recently discovered superconductivity of boron doped diamond 
within the framework of a phonon-mediated pairing mechanism. The role of the dopant, in substi- 
tutional position, is unconventional in that half of the coupling parameter A originates in strongly 
localized defect-related vibrational modes, yielding a very peaked Eliashberg a^F{ijj) function. The 
electron-phonon coupling potential is found to be extremely large and Tc is limited by the low value 
of the density of states at the Fermi level. 
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In a recent paper [1] , the superconductivity of boron 
doped (B-doped) diamond was evidenced experimentally 
with a transition temperature Tc of 4 K. This work fol- 
lows the discovery of superconductivity in doped silicon 
clathrates [2,3] (Tc=8 K), a cage-like silicon material 
which crystallizes in the same sp^ environment as the di- 
amond phase. Further, the superconductivity of carbon- 
based clathrates was predicted on the basis of ab initio 
simulations [3,4]. 

Even though the reported temperatures are rather low, 
the superconducting transition of column IV semicon- 
ductors is of much interest as it concerns very common 
materials. Further, in the case of the carbon clathrates, 
it has been predicted that the electron-phonon potential 
Vep=A/N(EF) (where N(Ei?) is the density of states at 
the Fermi level), may be extremely large [3,4]. Finally, 
the superconductivity of hole-doped strongly bonded co- 
valent systems is reminiscent of that of the MgB2 com- 
pound [5] . 

As compared to the cage-like clathrate phase, diamond 
is difhcult to dope as the network is much denser. This 
explains why the superconductivity of doped diamond 
was never observed, as only quite recently a doping in the 
limit of a few percent could be achieved. As described in 
Ref. [1], it is only for a boron concentration c > 10^^ cm~^ 
that a degenerate semiconductor or metallic character 
can be observed. 

In two recent work [6,7], a first theoretical study of 
the superconductivity of B-doped diamond within an el- 
egant virtual crystal approximation was proposed. In 
this approach, the crystal is made of "averaged carbon- 
bore atoms" (nuclear charge Z=(l-x)Zc + xZb) within 
the standard 2 atoms/cell structure of diamond, and B 
atoms are not considered explicitly. As such, the super- 
conductivity of B-doped diamond was analyzed in terms 
of Bloch-like delocalized electron and phonon "average" 
states. 

In this work, we study by means of ab initio calcula- 
tions the electron-phonon coupling in B-dopcd diamond. 



We consider explicitly a large supercell with one B atom 
in substitution. The metallic character of the doped sys- 
tem is evidenced and the deviations from a simple rigid 
band model are discussed. We show further that the 
boron atoms yield very localized vibrational modes which 
contribute to half of the electron-phonon coupling pa- 
rameter A, in large contrast with the picture conveyed 
by the virtual crystal approximation. The value of A is 
calculated ah initio^ yielding an unusually large electron- 
phonon coupling potential Yep or deformation potential 
D. 

Our calculations are performed within the local density 
approximation to density functional theory [8] (DFT) 
and a pseudo-potential [9] plane-wave approach. A 
55 Ryd energy cut-off and a (4x4x4) Monkhorst-Pack [10] 
sampling of the Brillouin zone (BZ) showed good con- 
vergence for structural relaxations. The phonon modes 
and the electron-phonon (e-ph) coupling matrix ele- 
ments were obtained within the framework of perturba- 
tive DFT [11]. Due to the computational cost, phonons 
were calculated at the q = F point only. This is equiva- 
lent to a (3x3x3) sampling of the undoped diamond BZ. 

For the calculation of the average e-ph matrix elements 
gqi, over the Fermi surface, required to evaluate A: 

A = N{Ef)Vep = 2N{Ef)Y, < > (1) 



< M > = 



(Pk 



BZ 



^ \gq^{knn')\ 



N{EfY 
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a much larger (8x8x8) k-point sampling was used. In 
Eqn. 2, the origin of the electronic energies (efc„) is taken 
to be the Fermi energy E^?. 

Our system consists of a diamond (3x3x3) supercell 
(54 atoms) with one B atom in substitution (noted BC53 
in what follows). We are thus in a 1.85% doping limit, 
close to the 2.5 ±0.5% studied experimentally. Upon 
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structural relaxation, the B-C bond length is 1.57 A as 
compared to the 1.54 A for the C-C bonds in undoped 
diamond (LDA values). The distance between second 
and third carbon neighbors is reduced to 1.52 A in order 
to accommodate the larger B-C bonds. 

We first explore the electronic properties and compare 
the band structure of BC53 with that of diamond (see 
Fig. 1). As expected, we are clearly in a metallic regime 
with the Fermi level located at ^ 0.54 eV below the top of 
the valence bands. Even though Fig. 1 seems to confirm 
the model of rigid band doping around E^, clear differ- 
ences exist in particular at the zone boundaries where 
degeneracies arc left due to the reduction of the peri- 
odicity under doping. We note that in most directions, 
the band structure of pure and B-doped diamond start 
to significantly depart from each other precisely below 
Ef [12]. 

We also calculate the total density of electronic states 

(eDOS) projected onto the atomic orbitals. Even though 
the eDOS projected onto the B orbitals (eDOSs) is 
small, it is instructive to plot (Fig. lb) the ratio 
NaXcDOSs/eDOS, where Na=54 is the number of atoms 
per unit cell. Clearly, close to E^, the relative weight of 
B orbitals into the wavefunctions is much larger than 
the weight on any given carbon atom. As expected, this 
ratio increases where the bands of the ideal and doped 
diamond differ most significantly. An analysis of the pro- 
jection of the eigenstates onto the atomic orbital basis 
shows that the leading coefficients are the {px,Py,Pz) B 
orbitals [13]. Such results indicate that even though de- 
generate with the diamond Bloch states, the wavefunc- 
tions of interest for e-ph coupling exhibit some degree of 
localization around the B atom. 

We now study the phonon states at F. The related 
DOS (pDOS) is shown in Fig. 2b and compared in Fig. 
2a to the one of undoped diamond [14]. As compared to 
the virtual-crystal approximation results, we do not ob- 
serve the very large softening of the zone-center phonon 
modes predicted in Ref. [6]. However, the splitting of 
the highest diamond mode (Fig. 2a) clearly results in 
transferring spectral weight to lower energy. This can 
be interpreted on the average as a softening of the opti- 
cal frequencies. In particular, zone-center modes emerge 
around 1200 cm^l (Fig. 2b) which are consistent with 
the appearance of new broad peaks in Raman experi- 
ments around 1220 cm~^ [1,15]. 

More insight can be gained by calculating the pDOS 
projected onto the B atom and its 4 C neighbors (thick 
line). The most salient feature is that the new modes 
(3-fold) at 1024 cm"! (vertical dotted hues in Fig. 2) 
project to more than 50% on these 5 atoms. The anal- 
ysis of the associated eigenvectors clearly indicates that 
these three modes correspond to the stretching of the 
B to neighboring C atom bonds. The displacement of 
the B atom is ^ 22% larger than that of C neighbors as 
expected from the difference in masses. Therefore, the 



dopant atoms induce vibrational modes which are very 
localized in space. This is a significant deviation from 
a virtual crystal approach where only extended phonon 
states can be evidenced. The influence of such modes on 
the e-ph coupling is now shown to be crucial. 

We follow equations 1-3 to calculate A. With our lim- 
ited q-point sampling [16], we find that X{q = F) is rather 
sensitive to the gaussian broadening used to mimic the S- 
functions in Eq. 2. However, such a variation is primarily 
related to the variation of the nesting factor n(q) (that 
can be obtained by setting gq,y{knn') to unity in Eq. 2). 
This factor n(q) is a measure of phase space availability 
for electron scattering by phonons. Indeed, one verifies 
that the quantity X{q)/n{q) (q=F) is much more stable 
and displays a plateau over a large energy broadening 
range. We thus define an averaged A by rewriting: 

A = y d\X{q) = j d\n{q) X{q)/n{q) 

~A(F)/n(F) J d^qn{q)=N{EFfX{r)/n{r) 

Such an approximation relies on the usual assumption 
that the e-ph matrix elements are varying much more 
smoothly as a function of q than the nesting factor. 

We find A=0.43 ±0.01. Such a value of A is much 
smaller than that of MgB2 (in the absence of anhar- 
monic effects [17], A(MgB2) ~ 1). Assuming now that 
the McMillan expansion [18] is valid in the limit of de- 
generate semiconductors: 



1.2kB 



exp 



-1.04(1 A) 



X - ^i*{i + omx) 



and with a logarithmic averaged phonon frequency LUiog 
of 1023 cm^'^, the screening parameter fi* needed to re- 
produce the 4 K experimental value of Tc is found to 
be of the order of 0.13-0.14 [19]. This is a very standard 
value, indicating that phonon-mediated pairing mecha- 
nisms can adequately account for the superconducting 
transition in impurity-doped diamond. 

With N(Ef)=20.8 states/spin/Ry/Cell, the e-ph po- 
tential Vep is 280 meV ± 10 meV. Such a value is 
comparable to what was found in the Li-doped car- 
bon clathrates {Vep ~ 250 meV [4]) and much larger 
than the 60-70 meV found for the fuUerides [20]. This 
was analyzed in terms of the amount of sp^ charac- 
ter of the C-C bonds [3]. For sake of comparison 
with Refs. [6,7], we calculate the deformation poten- 
tial D such that A=N2DVMa;2, where N2 is the eDOS 
in eV/state/spin/"2 atoms cell". With N2=0.057 and 
a;=1023 cm-\ one finds D = 27 eV/l-\ larger than 
the D ~ 21 eV/A~^ value found in the virtual crystal 
approximation [6,7]. This value is more than twice as 
large as the one calculated for MgB2, indicating that it 
is primarily the low density of states at the Fermi level, 
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related to the 3D nature of diamond, that is responsible 
for the low value of A and Tc- 

Better q-point sampling and the influence of anhar- 
monicity [6,17] will certainly change the value of A. Be- 
yond the numerical estimates, it is important to look at 
the spectral decomposition of A, that is the Eliashberg 
function a^F{uj) (Fig. 2c): 

a^F{ui) = ^^U}q^Xq^6{uJ - Wq^) (4) 

By comparison with the pDOS, it appears that a signifi- 
cant contribution to A comes from the B-related phonon 
modes. As a matter of fact, these 3 modes contribute to 
about 50% of the coupling [21]. In particular, the result- 
ing Eliashberg function is very peaked at the impurity vi- 
brational modes energy, displaying a (5-like shape. There- 
fore, not only the density of states N(Ef), but also the 
strength of the e-ph coupling, can be taylored by chang- 
ing the dopant concentration and/or chemical type. This 
mechanism offers thus an additional flexibility in trying 
to increase T^. The situation is very different from the 
one of fullerides or clathrates where the low frequency 
phonon modes associated with the dopant hardly con- 
tribute to A. 

To strengthen that point, we look at the superconduc- 
tivity of diamond within a rigid band model [22]. We 
plot in Fig. 3 the evolution of A as a function of the 
position of the Fermi level. For E/ located ~ 0.5-0.6 
eV below the top of the valence bands, we find a A of 
~ 0.24, that is roughly 55% of what we find in the BC53 
case. Even though very qualitative, this computer exper- 
iment substantiates the conclusion that B-related modes 
contribute significantly to the e-ph coupling. Further, 
Fig. 3 clearly suggests that larger A could be achieved 
under larger doping conditions due to the quadratic rise 
of N(Ei?) as a function ofEp- 

It is certainly accidental that just half of the coupling 
originates in B-modes as the weight of the eigenstates on 
the dopant should scale as y/Na under varying concen- 
tration (neglecting the electronic localization effect de- 
scribed above). In addition, the simple structural model 
studied here does not account for the effect of random- 
ness, nor the possibility for B atoms to cluster or adopt 
other configurations. It is clear however that the role 
of the dopant cannot be interpreted solely as provid- 
ing charges to the network, nor as softening the phonon 
modes by depleting the bonding states at the top of the 
valence bands. In particular, such results reveal that any 
attempt to explain the superconductivity of doped dia- 
mond within a rigid band model or the virtual-crystal ap- 
proximation should be taken with care. The importance 
of the impurity related vibrational modes is expected to 
be very general and should apply to any semiconduc- 
tor/insulator doped by strongly bound atoms. 

In conclusion, we have shown that the role of boron 



in doped diamond is unconventional in that a significant 
fraction of the coupling coefficient A originates in vibra- 
tional modes very localized on the defect centers. With 
a A of 0.43, much smaller than that of MgB2 ^ 1, Tc re- 
mains small. Even though the electron-phonon potential 
is extremely large, the 3D nature of the network reduces 
the density of states at the Fermi level. This invites to 
study the case of doped diamond surfaces where both 
the contraction of the reconstructed bonds and the 2D 
nature of the surface states may lead to much larger T^ 
[23]. 
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FIG. 1. (a) Details of the top of the valence bands along 
high-symmetry directions of the FCC BZ for undoped dia- 
mond in the (3x3x3) cell (dotted lines) and the doped BC53 
system (full lines). The horizontal line is the position of the 
Fermi level in the doped case. The two band structures have 
been aligned at the top of the valence bands, (b) Ratio of the 
density of states projected onto the B-orbitals to the total 
density of states per atom in the unit cell (this ratio has been 
set to zero in the band gap). 
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FIG. 2. Phonon density of states pDOS for (a) bare dia- 
mond with modes on an unshifted (3x3x3) grid (b) the boron 
doped BC53 cell with modes at P. In (c) the spectral function 
a^F(w) is represented for modes at P. Energies are in cm~^ 
and the y-axis magnitude is arbitrary. A 5 cm~^ broadening 
has been used. The vertical dotted lines indicate the energy 
position of the 3-fold B-related vibrational modes. 
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FIG. 3. Evolution of A for diamond, doped within a rigid 
band model (see text), as a function of the Fermi energy po- 
sition Ef- The energy reference has been set to the top of the 
valence bands. 
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